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Abstract 

We present a new method for obtaining information on the small x behavior 
of the structure function F2 outside the kinematic range of present acelerators 
from the mean inelasticity parameter in UHE neutrino-nucleon DIS interac- 
tions which could be measured in neutrino telescopes. 



1 Introduction 



A variety of models have been proposed in which astrophysical neutrinos exceed the 
expected atmospheric neutrino flux at energies in the TeV scale or higher The de- 
tectors presently in design or construction stages 0] look for Cerenkov light from the 
neutrino-induced muon in deep inelastic scattering (DIS) charged current (CC) interac- 
tions and take advantage of both the long muon range and the rise of the neutrino cross 
section to meet the neutrino detection challenge. 

Several alternatives have been proposed for observing high energy neutrinos and they 
are all based on the detection of the showers that are also produced in most of the neutrino 
interactions. These include the detection of particle shower fronts in the atmosphere in 
the horizontal direction (Horizontal Air Showers) |3j and coherent pulses from showers in 
different media, both acoustic || and in radio 0, |5[] waves. Shower detection is sensitive 
both to neutral current interactions and to all neutrino flavors. In DIS interactions with 
nuclei, showers of hadronic type are always initiated by the debris of the struck nucleons. 
For electron neutrino charged current interactions an electromagnetic shower is initiated 
at the lepton vertex in addition to the hadronic shower. Certainly conventional detectors 
in construction can also look for the Cerenkov light from the induced showers provided 
they are contained or sufficiently close to the instrumented volume. In fact it is likely that 
shower detection is the main way these detectors can observe ultra high energy neutrinos 
because of the earth's opacity and the atmospheric muon background. 

The high energy neutrino cross section is a crucial ingredient in the calculation of 
the event rate in high energy neutrino telescopes. The DIS cross section is given in 
terms of structure functions whose energy dependence (scaling violations) is calculated in 
perturbative QCD. In practice, structure functions are computed from available parton 
(quark and gluon) distribution functions which depend on x and Q 2 (see below). Parton 
densities are extracted from data using DGLAP evolution equations which effectively 
sum ln(Q 2 ) terms contained in the QCD perturbative expansion. 

The DGLAP approach can break down at low x because of potentially large ln(l/x) 
terms which also appear in the perturbative series. The theoretical behaviour of structure 
functions at small x has been discussed since the seventies and recently revived in the 
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context of the low x HERA data (see Ref. J7| for a review). 

The sumation of the ln(l/x) terms for intermediate Q 2 is performed with the BFKL 
formalism H which evolves parton densities unintegrated in transverse momentum. Within 
BFKL at LO and fixed a s , the structure function F 2 behaves as power like x~ x , with 
A ~ 0.53 for a s =0.2 (for a s (M 2 )=0.12, A ~ 0.32). The exponent A is related to the inter- 
cept of the Pomeron which governs the high-energy asymptotics of the total cross section. 
The consideration of NLO effects in BFKL is a matter of present discussion due to the 
large correction found (A becames negative). However, recently it has been found |J that 
the NLO value of A, improved by optimization of the renormalization scale, has a very 
weak dependence on the virtuality Q 2 , ranging the values A ~ 0.13 — 0.18 at Q 2 = 1 — 100 
GeV 2 . 

This result, x~ x , is obtained in DGLAP when the input structure function at fixed Qq 
has a singular behavior at low x but also when it is flat. In the latter case, the solution is 
the so call double logarithmic approximation (DLA) and effectively resums terms of type 
o; s ln(l/x) ln(Q 2 ) [PHI - After evolution over a suficiently large Q 2 interval, the DLA result 



mimics the BFKL behavior, x x ill 



At this point it should be stressed that the low x behavior of structure functions is 
fundamental in the high energy limit of the neutrino-nucleon DIS cross section. The 
calculation of the neutrino-nucleon cross section involves an integration over x which 
corresponds to the fraction of momentum carried by the struck parton (see below). For 
x~ x with A > 0, as the neutrino energy increases, the integral becomes dominated by 
the interaction with partons of lower x, while the Q 2 integral remains dominated by Q 2 
values up to the electroweak boson mass squared (see for example |L2[ and below). For 
Q 2 above the integrand behaves as Q~ A and quickly becomes irrelevant. 

On the experimental side there are however no F 2 measurements at very small x and 
large Q 2 . Consequently, in the process of the total cross section calculation, the parton 
densities extracted from present data must be carefully extrapolated to the region in the 
x, Q 2 plane without experimenta support. The dominant kinematical region calls for the 
connection between the low x behavior at intermediate Q 2 (the BFKL region) and the 
high Q 2 region but at moderate low x obtained in fits to experimental data (the DGLAP 
region) JIE| . 
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The consideration of QCD effects in the neutrino DIS cross section for neutrino en- 
ergies well above the electroweak boson threshold was done at the end of the seventies 
13|, |T4| motivated by the DUMAND project. Previously the neutrino-nucleon DIS cross 



section had been studied in detail for lower neutrino energies in an attempt to disentangle 
asymptotic freedom effects and, in particular, scaling violations (see |15) for a review). 
The study of the uncertainty in the cross section calculation, which is important in the 



context of high energy neutrino telescopes, has been recently addressed |17|, [18], [19| . The 
uncertainty due to the extrapolation to high Q 2 is not expected to be large because the 
low x shape of the parton distribution functions at low Q 2 is narrowly constrained by 
HERA while the further evolution to high Q 2 values seems well defined by DGLAP QCD 
evolution. At the highest energies uncertainties within 20% JlT], 40% |TI| and a factor 



2 ±:L |18| are typically reported. However we think that the consideration of conventional 
DGLAP evolution to small x values and high Q 2 , as for example x = 10~ 9 and Q 2 = 10 4 
GeV 2 , should at least be questioned in view of absence of data. 

Conversely if any neutrino interactions measurements could be made leading to infor- 
mation on the cross section, such information would be of great value to constrain the 
theoretical predictions at low x and high Q 2 . Of particular interest is the y distribution 
which can be shown very sensitive to the value of A (see below), y being the fraction of 
the neutrino energy which flows to the hadronic part of the interaction in the laboratory 
frame. 

This distribution is incidentally very important for neutrino detection, particularly for 
all the alternatives in neutrino detection which rely on detecting the high energy showers 
which are always produced at the hadronic vertex, whatever the neutrino flavor and for 
both charged current (CC) and neutral current (NC) interactions. 

The article is organized as follows: Firstly we collect the most relevant formulae used 
in the calculation of the neutrino nucleon cross section. Then, we present the result 
for the total cross section and the average inelasticity (y) computed from two sets of 
parton distribution functions. We explain the characteristics of the results for neutrinos, 
antineutrinos in CC and NC interactions in terms of the parton distribution dependences, 
stressing the connection between these quantities and the low x limit of the structure 
functions at the highest energies. 
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In the assumption that the low x behavior of F2 at large Q 2 can be described by x~ x , 
we present a very simple analytical relation between A and the average inelasticity (y) 
which does not depend of any other parameter. We suggest that A could be determined 
from y-measurements in the events which are expected in high energy neutrino detec- 
tors. Although ^-measurements could be thought to be rather speculative at this early 
stage neutrino astronomy is at present, we point out a method to measure y in neutrino 
telescopes. 

2 The neutrino-nucleon DIS cross section at UHE 

In terms of structure functions the charged-current (CC) neutrino-nucleon DIS differ- 
ential cross section is given by: 

iSr = &) 2ME - (sfe) 2 ^ ~ ^ ± « 

where y± — 1 ± (1 — y) 2 , M is the nucleon mass, E v the neutrino energy in the lab frame, 
Q 2 = 2ME u xy and we have neglected terms suppressed by powers of M 2 /Q 2 . 

In the QCD improved parton model, the structure functions Fi {i = 2, 3, L) are cal- 
culated in terms of quark and gluon distribution functions. At leading order (LO) ap- 
proximation in perturbative QCD F 2 {xF 3 ) is simply related to the sum (difference) of 
parton densities. At next-to-leading order (NLO) further integrals are involved relating 
at order a s (Q 2 ) the parton densities with the structure function. The contribution due to 
F L) which is also proportional to a s (Q 2 ) becomes small at large Q 2 and like the M 2 /Q 2 
terms are neglected in this work. 

We present, for simplicity, the expressions which relate the cross sections and the par- 
ton densities at the LO approximation. Within the above approximations the differential 
cross section for the charged current (CC) interaction of neutrinos on isoscalar target 
takes the form: 



d 2 a uN Gl M^ME, 



dxdy 7T (M^ + Q 2 )' 

d 2 a pN Gl MhrME v 



dxdy 7T (M^ + Q 2 )' 



A(x,Q 2 ) + B(x,Q 2 )(l-yf 
A(x,Q 2 )+B(x,Q 2 )(l-y) 2 (2) 



where 

A{x, Q 2 ) = x(u + d + 2s + 2b) , B(x, Q 2 ) = x(u + d + 2c + 2t) 

The antiquark combinations A and B are obtained from A and B replacing each quark 
by the corresponding antiquark of the same flavor. 

The total cross section is calculated by integration of Eq. (2). We have used two 
representative sets of parton densities^: the LO set from Ref. |^TJ (MRST98) and the 
NLO set in the DIS factorization scheme from Ref. |22|] (GRV98)0. In the integration we 
extrapolate the quark distribution functions in x below the low x limit given by the authors 
(x = 10~ 5 for MRST98 set and x = 10~ 9 for the GRV98 set) using the x slope at the 
lowest x value of the parametrization for each value of Q 2 . This simple phenomenological 
extrapolation agrees with the more elaborated prescriptions based on perturbative QCD 
as for example the double-logarithmic-approximation which explain HERA data (see [JEJ 
and also Fig. lb). We choose this ad hoc extrapolation for simplicity because our purpose 
is just to show the sensitivity of the high energy neutrino cross section to the low x parton 
behavior. 

The total charged and neutral current interaction cross sections^ for both neutrinos 
and anti-neutrinos calculated with MRST98 partons are shown in Fig. la as a function 
of the neutrino energy. At low energy one can observe the linear rise of the cross section 
with the neutrino energy (E„). At high energy it would follow a logarithmic rise \n(E u ) 
were it not for the increase of the sea quark densities at small x (~ x~ x ) which imply a 
high energy cross section that rises as (E v ) x ln(E u ) for E v above 1 PeV. 

We have also calculated the cross sections with GRV98 partons. For neutrino energies 
considered in this work, i.e. below 10 13 GeV, the results are insensitive to the extrapolation 
below x = 10~ 9 . The cross sections from MRST98 and GRV98 parton distribution can 
be compared in Fig. lb. 



1 The consideration of nuclear effects in parton distributions which could be relevant for heavy nuclei 
will be presented elsewhere p(| 

2 Rigurously, NLO corrections modify Eq. (2) in the DIS scheme by terms proportional to a s coming 
from XF3. We neglect them in this calculation 

3 For brevity, in this work we do not present the explicit expressions used in the calculation of the NC 
cross sections 
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3 The mean inelasticity at UHE 



In this work we are interested in the fraction of the neutrino energy which flows to the 
hadronic part of the interaction in the laboratory frame, y, which is also called inelastic- 
ity. For a given muon neutrino high energy flux, this parameter fixes the relative rates of 
the two main types of detections, using muons in charged current interactions and using 
the showers produced in the interactions. It is also reponsible for the relative sizes of 
the electromagnetic and hadronic showers induced in charged current electron neutrino 
interactions. This quantity is necessary for any attempt at extracting the neutrino en- 
ergy in high energy neutrino telescopes whether by detecting the shower produced in the 
interaction or by detecting the muon produced by muon neutrino CC interactions from 
the detected hadronic shower or muon. 

The average value of y can be obtained by integration of the differential cross section: 

(y) = -[ 1 dyy^ (3) 
a Jo ay 

The energy dependence of (y) has been studied in the past to test asymptotic freedom in 
strong interactions (see |]TSJ for a review). The effect of the boson propagator on the y 



distribution and the relevance of the small x region were pointed out in Refs. [13, 14 



In the present work we have calculated (y) using modern parton densities which better 
describe the low x kinematic region. The results using MRST98 partons are shown in 
Fig. 2a. 

Let us explain the general features observed in Fig. 2a from the cross section formulas 
given in Eq. (2). At low energy, the value of (y) for v is larger than for v because the most 
important contribution comes from the valence quarks which are suppressed by (1 — y) 2 
in the anti-neutrino cross section formula (see Eq. (2)). If there were not sea quarks, the 
y distributions given by Eq. (2) would behave as (1 — y) 2 (would be flat) for antineutrinos 
(neutrino) and correspondingly (y) = 0.25 (0.5). In the case of low energy z/, the increase 
of (y) with energy is due to the rise of the sea quark distribution functions. For neutrinos, 
this effect is less apparent. 

The depletion of (y) for energies above 1000 GeV (see Fig. 2a) is due to the W 
propagator appearing in Eq. (2) which acts as a cutoff in the integration restricting the 
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values of Q 2 to around M^, i.e. xy ~ My V -/{2ME v ). Thus, the da/dy distribution is 
shifted to lower values of y. 

At the highest energies, because of the sea quark dominance, both neutrino and an- 
tineutrino cross sections become almost equal (see Fig. 2a) and the small difference 
between CC and NC interactions is due to the characteristic couplings of the electroweak 
interaction. 

Concerning the value of (y), Fig. 2a shows that it typically becomes stable at values 
slightly above 0.2. It is remarkable that in the high energy limit (y) should be constant 
while the average value of the Bjorken x variable, (x), decreases strongly (see Fig. 4b). 
The low x rise of the sea quark densities (predicted by perturbative QCD) shifts the 
weight of the cross section to small x values compensating the y distribution. 

The energy dependence of (y) at high energy is related to the behavior of the slope 
assumed constant in the extrapolation of the MRST98 quarks at low x. If A were negative 
or zero, (y) would send to zero in this limit. 

For example, if in the cross section calculation one uses GRV98 partons, (y) is not 
constant and we find that it decreases with increasing energy (see Fig. 2b) which reflects 
the depletion of the effective slope A of the GRV98 partons with increasing Q 2 . 

In the following we show that the observed correlation between (y) at high energy and 
A at low x can be put in analytic form using very simple approximations. 

4 Low x physics and high energy neutrino telescopes 

Let us assume that F 2 at low x and high Q 2 is given by the power-like expression 
F 2 = A(Q 2 ) x~ x , where A is assumed constant (or smoothly dependent on Q 2 as predicted 
by ||) Then the CC neutrino-nucleon differential cross section (in y and Q 2 ) takes the 
form: 

dydQ 2 -\A,)\M 2 w + Q 2 ) A{Q) { Q 2 ) y V (4) 

where in Eq. (1) we have neglected the contribution from F L and xF 3 structure functions, 
expected to be negligable at low x and high Q 2 . 

In the calculation of (y) by Eq. (3) the Q 2 integral of Eq. (4) in numerator and 
denominator cancels at high energy, i.e. provided that 2ME U is much larger than M^. 



Then one has: 

= &dyy% J^dy [1 + (1 - yf] 

W Iody% Jo 1 dy [1 + (l - yf] y^ {) 

which can be easily integrated to get the simple analytical relation: 

_ A 3 + 5A 2 + 8A 
{V) ~ A 3 + 6A 2 + 13A + 12 U 

Analogously, one can proceed as above in the calculation of the average Q 2 which is 
expected (see above) to be around at energies sufficiently above the W propagator 
threshold. One finds: 

To check Eq. (6), let us find out the value of A which corresponds to (y) = 0.23, which 
is the result from the explicit integration of the CC neutrino-nucleon differential cross 
section da/dy calculated with MRST98 parton densities at E v = 10 10 GeV (see Fig. 2a). 
Looking at Fig. 3 one finds A = 0.42. We have checked that the MRST98 sea partons at 
x = 10 -5 and Q 2 = 10 4 GeV 2 have this slope which confirms the validity of Eq. (6). 

Furthermore , substituting A = 0.42 in Eq. (7), one obtains (Q 2 ) = 0.9 10 4 GeV 2 
which agrees with the asymptotic value reached at high energy from explicit numerical 
integration of the diferential cross section using the MRST98 partons (see Fig. 4a). 

Assuming that the mean inelasticity (y) could be experimentally determined, Eq. (6) 
can be inverted to obtain directly the F 2 slope parameter A at small x for Q 2 around M^, 
which is outside the kinematic limits of present accelerators. 

We have to stress that the measurement of (y) will not be easy to obtain because 
it would require the determination of the final state energy corresponding to both, the 
nuclear cascade and the lepton. If electron neutrinos are detected, the showers that come 
out of charged current interactions will have a different character depending on y because 
the neutrino induced electromagnetic (hadronic) shower arising at the leptonic (hadronic) 
vertex carries a fraction 1 — y (y) of the neutrino energy. In particular it has been recently 
suggested that it may be possible to measure y by using the radio technique ||. This 
tecnique relies of the detection of coherent radio pulses from the excess charge in the 
induced showers in a dense medium such as ice. As the radiation is coherent the angular 
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structure of the distribution is sensitive to the longitudinal development of the shower. 
The sensitivity to y arises from the elongation of the showers that develop at the electron 
vertex because of the Landau-Pomerancuck-Migdal (LPM) |24| effect ||. 

For the case of muon neutrinos in charged current interactions, both the muon and the 
nuclear shower in the final state would have to be detected. As the muon energy loss is 
proportional to the muon energy above 1 TeV, this may not be out of question in a future 



large conventional neutrino telescope such as IceCube |23[ although many difficulties can 
be forseen. In any case the ratios of rates detected with muons to those detected with 
showers should be dependent on y for a given high energy neutrino flux. 

5 Conclusions 

The measurement of (y) at high energies in high energy neutrino interactions would 
give direct information on the behavior of the sea quarks distribution functions in the 
nucleon in a kinematic range unreachable to the most powerful accelerators, independently 
of the neutrino flux. We have presented a simple analytical parametrization which relates 
(y) with A, the x slope of F 2 at low x and high Q 2 . 

We have also pointed out that neutrino telescopes are sensitive to the value of (y) 
and the possibility to measure y from the interference pattern of the radio-signals of the 
showers produced in the neutrino interaction. 

Some other interesting information can be obtained from the measurement of y. For 
example, it has been proposed that supersymmetric effects, as i?-parity violation could 



increase {y) with energy significantly [ 25| . Also the effect of extra dimensions in neutrino- 



nucleon interactions [21) ] should be manifested through the modification of (y) p7| . 
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Figure captions 



Figure 1. The total neutrino-nucleon deep inelastic cross section as a function of 
neutrino energy in lab frame for a) NC and CC neutrino and antineutrino interactions 
with MRST98 partons and b) CC neutrino interactions with MRST98 and GRV98 parton 
densities. 

Figure 2. The average inelasticity y in deep inelastic neutrino-nucleon interaction 
as a function of laboratory neutrino energy for a) NC and CC neutrino and antineutrino 
interactions with MRST98 partons and b) CC neutrino interaction with MRST98 and 
GRV98 parton densities. 

Figure 3. The parameter A (F 2 ~ x~ x ) at low x and high Q 2 ~ as a function of 
the average inelasticity in CC uN deep inelastic interaction. 

Figure 4. a) The average Q 2 in deep inelastic CC neutrino-nucleon interaction as 
a function of the laboratory neutrino energy, b) The average x in deep inelastic CC 
neutrino-nucleon interaction as a function of the laboratory neutrino energy. 
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Fig. la 
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Fig. lb 
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